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Abstract. This is a brief review of the standard model of cosmology. We first introduce the FRW 
models and their fiat solutions for energy fluids playing an important role in the dynamics at different 
epochs. We then introduce different cosmological lengths and some of their applications. The later 
part is dedicated to the physical processes and concepts necessary to understand the early and very 
early Universe and observations of it. 
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The purpose of the present review is to provide the reader with an outline of modern cos- 
mology. This science is passing through a revolutionary era, mainly because recent high 
precision observations have severely constrained theoretical speculations, and opened 
^ new windows to the cosmos. Tracking the recent history, in the late 1940's George 

Gamow [1] predicted that the Universe should had begun from a very dense state, char- 
acterized by a huge density at very high temperatures, a scenario dubbed the Big Bang, 
that was conjectured by George Lemaitre in the early 30's. This scenario predicts that 
matter and light were at very high energetic states, and both components behaved as a ra- 
diation fluid, in thermal equilibrium described as Planck blackbody (Bose-Einstein and 
the related Fermi-Dirac distributions). This initial state remains today imprinted in the 
Cosmic Microwave Background Radiation (CMBR). Gamow's scenario predicted this 
primeval radiation would be measured at a temperature of only few Kelvin's degrees; 
since the expansion of the Universe cools down any density component. Robert Dicke, 
and others, begun the race to discover this radiation coming from the cosmos, and in 
1965 A. A. Penzias and R.W. Wilson of the Bell telephone laboratories discovered, by 
chance, this radiation form, confirming the general Big Bang scenario; see the original 
references published in Ref. [2]. 

During this elapsed period the Big Bang scenario was generally accepted. However, 
some key questions remained open, for instance, whether or not this radiation was of 
Planckian nature to entirely confirm that the Universe was in thermal equilibrium at the 
very beginning of time. In 1990 a modern version of the Penzias-Wilson experiment 
was carried out. This experiment, lead by G. Smoot and J. Mather with the Cosmic 
Background Explorer (COBE) satellite, started a new high precision experimental era 
in cosmology [3]. The COBE team for the first time revealed that the Universe was 



in equilibrium (the radiation had the Planck form to high precision), and was almost 
homogeneous and isotropic, but not completely [4]. The tiny (10~ 5 ) anisotropics found 
by COBE -also imprinted in the matter distribution- were lately the responsible for the 
formation of stars, galaxies, clusters, and all large scale structures of our Universe. 

The origin of these tiny anisotropics is presumably in quantum fluctuations of funda- 
mental fields of nature which were present in the very early Universe. Modern quantum 
field theories, together with cosmological models, help to understand how this small 
fluctuations evolved to become of cosmic scales, enabling COBE to detect them. This 
satellite and other more recent cosmological probes, such as BOOMERANG, MAX- 
IMA, WMAP and now PLANCK satellite that measured CMBR and the Two degree 
Field (2dF) galaxy survey and the Sloan Digital Sky Survey (SDSS), among others, not 
only confirmed with a great accuracy some of the theoretical predictions of the standard 
big bang cosmological model, but also opened the possibility to test theories and sce- 
narios applicable in very early Universe, such as inflation, or in present times, such as 
quintessence. In this way, cosmology that used to be a purely theoretical and often spec- 
ulative science, is today subject to high precision tests in light of these new observations 
[5,6]. 

ON THE STANDARD MODEL OF COSMOLOGY 

We begin our study by reviewing some aspects of the standard lore of physical and the- 
oretical cosmology. In doing that we consider the Friedmann-Robertson- Walker (FRW) 
model in Einstein's general relativity (GR) theory. We shall make use of "natural" units 
h = c = kg = 1 and our geometrical sign conventions are as in Ref. [7]. 



FRW models 

The cosmological principle states that the Universe is both spatially homogeneous and 
isotropic on the large scale, which was originally assumed but now observed to be valid 
for the very large large scale of the Universe. This homogeneous and isotropic space- 
time symmetry was originally studied by Friedmann, Robertson, and Walker (FRW), see 
Refs. [8]. The symmetry is encoded in and defines the unique form of the line element: 



ds 2 = -dt 2 + a 2 (t) 



dr 2 

= + r 2 {dQ 2 + sin 2 d<p 2 ) 

1 — kr z 



(1) 



where t is the cosmic time, r — 9 — (j) are polar coordinates, which can be adjusted so 
that the constant curvature takes the values k = 0, +1, or — 1 for a flat, closed, or open 
space, respectively. a(t) is the unknown potential of the metric that encodes the size at 
large scales, more formally is the scale factor of the Universe. 

To then solve this line element for the scale factor of the Universe a(t) and perturba- 
tions one simply needs a valid theory of gravity. The standard model of cosmology is 



based on GR, which can be derived from the Einstein-Hilbert Lagrangian 

(2) 
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where R is the Ricci scalar, G the Newton constant, and g = \g^ v | the determinant of the 
metric tensor. By performing the metric variation of this equation, one obtains the well 
known Einstein's field equations 

tf M v - ^ R 8nv = ^GT^ , (3) 

where T^ v = — j== ^ v 8 is the stress energy-momentum tensor, which is a symmetric 
tensor. Thus, Eq. (3) is a collection of ten coupled partial differential equations. How- 
ever, the theory is diffeomorphism invariant, and one adds to them a gauge condition, 
implying in general four extra equations to Eq. (3) that reduce the degrees of freedom. 

The beautifully symmetric FRW solutions to the Einstein Eqs. (3) represent a cor- 
nerstone in the development of modern cosmology, since with them it is possible to 
understand the expansion of Universe as was realized in the 1920's through Hubble's 
law of expansion [9]. With this FRW metric, the GR cosmological field equations are, 

(4) 
\a j j a" 

and 

h. AttCt 

(5) 
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where H is the Hubble parameter, p and p are the density and pressure of the perfect 

fluid considered, that is, T^ v = pu^u v + p(u^u v + g^. v ), where u^ = 8^ is the fluid's 

four velocity in comoving coordinates. Dots stand for cosmic time derivatives. 

Assuming energy-momentum tensor conservation, r„ v ;v = 0, is valid, one obtains the 

continuity equation, 

p + 3H(p+p)=0 . (6) 

The above system of equations (4), (5), and (6) implies three unknown variables (a, 
p, p) for three equations, but the equations are not all linearly independent, just two of 
them. Thus, an extra assumption has to be made to close the system. The answer should 
come from the micro-physics of the fluids considered. For the moment let us assume a 
barotropic equation of state that is characteristic for different cosmic fluids, w = const., 
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for radiation or relativistic matter 



P _ _ I for dust m 

p ~ W ~ | 1 for stiff fluid ( ' ' 



— 1 for cosmological constant 



to integrate Eq. (6), yielding 



M w 



fl 3(i+w) p \ a J 



where M w is the integration constant and is different dimensioned by considering differ- 
ent w— fluids. With this equation the system is closed and can be solved once the initial 
conditions are known. 

In addition to the above fluids, one can include an explicit cosmological constant (A) 
in Eq. (4) and arrange it in the following form: 

Q = n*+n M +n A = i + -^-2 ( 9 ) 

a L H L 

with Q.r being the radiation component dominated by the CMB and small (5 x 10~ 5 at 
present), Q.m = \Jr an d £2 a = thi ■ The parameter Q. is called the density parameter, 
which is composed of a matter part and a cosmological constant term. Thus, we see that 
the different values of the density parameters (Gl m , Cl\) will impose different values for 
the curvature term. If Q. > 1, it turns out k is greater than zero, signifying a Universe 
with a positive curvature (closed Universe). If Q. < 1, then k < 0, this corresponds to a 
negative curvature (open Universe). A critical value is obtained obviously when Q = 1, 
then the spatial curvature is null, k = 0. The value of the energy density for which 
Q. = [p +A/(SkG)]/p c = 1 holds is known as the critical density, p c = 3H 2 /8nG. 
The last term in Eq. (9) can be defined as £1% = — 2772 an d thus the Friedman equation 
becomes a constriction for the density parameters J^-Q,- = 1, and this expression holds 
at any time. It is useful to define an alternative measure of the expansion of the Universe 
through the redshift (z), 1 + Z = ao/a(t), where ao is the scale factor at present and is the 
set to unity by convention. Today zo = and in the early Universe the redshift grows. In 

terms of the redshift the density parameters are, from Eq. (8), Cli = Qr { (1 +z) 3 ^ l+w '', 
quantities with a subindex or superindex "0" are evaluated at the present time. 

Let us very briefly recall which w— values are needed to describe the different epochs 
of the Universe's evolution. The assumption that w — 1/3 is valid for a "fluid" of 
radiation and/or of ultra-relativistic matter (T ^$> m,m being its rest mass). This epoch is 
of importance at the beginning of the hot big bang theory, where the material content of 
the Universe consisted of photons, neutrinos, electrons, and other massive particles with 
very high kinetic energy. After some Universe cooling, some massive particles decayed 
and others survived (protons, neutrons, electrons) whose masses eventually dominated 
over the radiation components (photon, neutrinos) at the equality epoch (p K \ — p m ) 
at Zeq ~ 3200 [10]. From this epoch and until recent efolds of expansion (zde ~ 1/2) 
the main matter component produced effectively no pressure on the expansion and, 
therefore, one can accept a model filled with dust, w = 0, to be representative for the 
energy content of the Universe in the interval 3200 < z < 1/2. The dust equation of state 
is then representative of inert, cold dark matter (CDM). Dark matter (DM) does not 
(significantly) emit light and therefore it is dark. Another possibility is that dark matter 
interacts weakly and is generically called WIMP (Weakly Interacting Massive Particle), 
being the neutralino the most popular WIMP candidate. Another popular dark matter 
candidate is the axion, a hypothetical particle postulated to explain the conservation of 
the CP symmetry in quantum chromodynamics (QCD). Back to the Universe evolution, 
from z ~ 1/2 and until now the Universe happens to be accelerating with an equation 
of state w« -1, due to some constant energy that induces a cosmological constant, 
A = SnGp =const. The cosmological constant is the generic factor of an inflationary 



solution, see the k = solution below, Eq. (12). The details of the expansion are 
still unknown and it is possible that the expansion is due to some new fundamental 
field (e.g. quintessence) that induces an effective A(t) ~const. One calls (M. Turner 
dubbed it) dark energy (DE) to this new element. Dark energy does not emit light 
nor any other particle, as so far known, it simply behaves as a (transparent) media 
that gravitates with an effective negative pressure. The physics behind dark energy or 
even the cosmological constant is unclear since theories of grand unification (or theories 
of everything, including gravity) generically predict a vacuum energy associated with 
fundamental fields, < 0|r^ v |0 >=< p > g^ v , that turns out to be very large. This can 
be seen by summing the zero-point energies of all normal modes of some field of mass 
m, to obtain < p >ps M 4 /(l6n 2 ), where M represents some cutoff in the integration, 
M ^> m. Then, assuming GR is valid up to the Planck (PI) scale, one should take 
M ps 1/y/SnG, which gives < p >= 10 71 GeV 4 . This term plays the role of an effective 
cosmological constant of A = SnG < p >~ Mp, ~ 10 38 GeV 2 which must be added to 
the Einstein equations (3) and yields an inflationary solution Eq. (12). However, since 
the cosmological constant seems to dominate the Universe dynamics nowadays, one has 
that 

A ps SKGpo = ml ~ l(T 83 GeV 2 . (10) 

which is very small compared with the above value derived on dimensional grounds. 
Thus, the cosmological constraint and theoretical expectations are rather dissimilar, by 
about 121 orders of magnitude! Even if one considers symmetries at lower energy scales 
the theoretical A is indeed smaller, but never as small as the cosmological constraint: 
Aqut ~ 10 21 GeV 2 , A su n) ~ 10~ 29 GeV 2 . This problem has been reviewed since 
decades ago [11, 12] and remains open. 

The ordinary differential equations system described above needs a set of initial, or 
alternatively boundary, conditions to be integrated. One has to assume a set of two 
initial values, say, (p(f*),a(f*)) = (p*,d*) at some (initial) time t*, in order to determine 
its evolution. The full analysis of it can be found in many textbooks [13, 7]. Here, in 
order to show some physical, early Universe consequences we assume k = 0, justified 
as follows: From Eqs. (4) and (8) one notes that the expansion rate, given by the Hubble 
parameter, is dominated by the density term as a(t) — >■ 0, since p ~ 1/a 3 ^ 14 ~ w > > k/a 2 
for w > — 1/3, that is, the flat solution is very well fitted at the very beginning of times. 
Therefore, assuming k = 0, Eq. (4) implies 

a(t) = [6KGM w (\+wf}^)(t-h)^) 

(fnGMi) 1 / 4 (t - u) [ / 2 for w=\ radiation 

(67tGM o y/ 3 (t-u) 2 / 3 for w = 0dust (H) 

(lAnGMx) 1 / 6 (t-u) 1 / 3 for w = 1 stiff fluid 

and 

a(t) = a*e * for w = — 1 cosmological constant (12) 

where the letters with a subindex "*" are integration constants, representing quantities 
evaluated at the beginning of times, t = t*. To obtain Eq. (12), the argument given 



right above to neglect k is not anymore valid, since here p = const.; that is, from the 
very beginning it must be warranted that H 2 rs ^pp* > k/a 2 ,, otherwise k cannot be 
ignored. Nevertheless if A is present, it will eventually dominate over the other decaying 
components, this is the so called cosmological no-hair theorem [14]. A general feature of 
all the above solutions is that they are expanding, at different Hubble rates, H = jrnr^i \ 
for Eqs. (11) and H =const. for Eq. (12). 

From Eq. (11) one can immediately see that at t = t*, a* = and from Eq. (8), p* = °°, 
that is, the solution has a singularity at that time, at the Universe's beginning; this initial 
cosmological singularity is precisely the big bang singularity. As the Universe evolves 
the Hubble parameter goes as H ~ l/t, i.e., the expansion rate decreases; whereas 
the matter-energy content acts as an expanding agent, cf. Eq. (4), it decelerates the 
expansion, however, asymptotically decreasing, cf. Eqs. (5) and (8). In that way, H 
represents an upper limit to the age of the Universe; for instance, H = It for w = 1/3 
and H = 3t/2 for w = 0, t being the Universe's age. 

The exponential expansion (12) possesses no singularity (at finite times), being the 
Hubble parameter a constant. A fundamental ingredient of this inflation is that the right 
hand side of Eq. (5) is positive, a > 0, and this is performed when p ± 3p < 0, that is, one 
does not have necessarily to impose the stronger condition w = — 1, but it suffices that 
w < —1/3, in order to have a moderate inflationary solution; for example, w = —2/3 it 
implies a = aj 2 , a mild power-law inflation. 

Since the scale factor evolves as a smooth function of time (most of the time!), one 
is able to use it as a variable, instead of time, in such as a way that d/dt = aHd/da. 
This change of variable helps to integrate the continuity equation for non-constant w(a) 
to obtain: 

p(a)=p e~ 3/[1+w(a)lda/fl . (13) 

If, for instance, one parametrizes dark energy through an analytic function of the scale 
factor, w(a), one immediately obtains its solution in terms of 

1 ^ (14) 



y/SnGp(a)/3 a ' 

In cosmology, typical times and distances are determined mainly by the Hubble 
parameter, and in practice measurements are often related to redshift, as measured 
from stars, gas, etc. It is then useful to express the Friedmann Eq. (4) in terms of the 
redshift. The standard model of cosmology considers a Universe filled baryons, photons, 
neutrinos, CDM, and a cosmological constant (A), and is termed ACDM for short. For 
this model one obtains: 

H 2 = H 2 ^f\\+z? {i+Wi \ (15) 

i 

where W[ is the equation of state parameter for each of the fluids considered. The present 
contribution of the main energy components can be fitted from different cosmological 

probes, obtaining ^ 0) = 0.046 ±0.002, Q^m = 0.23 ±0.01, and d^ = 0.73 ±0.02, 
together with a Hubble constant of 70.4+}'4 km/s/Mpc [10]. Around 96% of mater- 
energy the Universe is composed of dark components! 



In general, if dark energy is a function of the redshift, from Eq. (13) one can generalize 
the above equation to: 

ff(z) 2 /Hj = nL 0) (H-z) 3 + 4 0) (l+z) 4 + «f(l+z) 2 + nSf(z), (16) 

where m stands for dark matter and baryons, 7 for photons, and 



f (z) = exp 



1+z 



(17) 



Eq. (14) gives the age of the Universe in terms of the redshift, Hq, and the density 
parameters: 

f°° dz 
to=H [ (18) 

Jo (l+z)H{z) 

When combining different cosmological probes one obtains for the ACDM model an 
ageof? = 13.75 ±0.11 Gyr [10]. 



Cosmic distances and their measurements 

It is useful to write the FRW metric, Eq. (1), in terms of a new distance coordinate (%) 

ds 2 = -dt 2 + a 2 (t) [dx 2 + fl (X)(d6 2 + sin 2 6 # 2 )] , (19) 

where 



sin*, 


* = +l, 


X, 


* = 0, 


sinh^ , 


k = -l. 



fk(x) = < X, * = 0, (20) 



Now we proceed to define some cosmic distances necessary to understand the cosmic 
physics. 

Causal horizon. The region of space that can be connected to some other region by 
causal physical processes, at most through the propagation of light, implies ds 2 = 0. For 
the FRW Eq. (1) or (19), in spherical coordinates with 0,0 =const. implies that [15, 13]: 

fXH [to dt' 1 f z dz' 

Xh= dx= -j^ = —^ \ t^tt < 21 ) 

Jo it ait') ciqHo Jo E(z') 

this is the so-called comoving distance, being the distance between two points in the 
Universe in which the expansion is factored out. Now, the causal or particle horizon, dn 
is given by: 

fXn ft dt' 

d H (t) = a(t) dx = a(t) —. (22) 

Jo Ju a \t ) 

In order to analyze the whole horizon evolution, from nowadays (to) to the Planck time 
(tpi), we have to consider all the Universe stages, but for brevity we shall not include the 



current accelerated expansion. We firstly compute the horizon for the matter dominated 
era t eq , < t < to and secondly for the radiation era t < t eq ,, because they are differently 
determined by Eq. (11), where we set t* = for convenience. For the matter epoch one 
has a{t) = ao(t/to) 2 ' 3 , then the first equation above gives Xh = ~(/o > fr° m the 
second equation one obtains the horizon djj{t) = 3t = 2H ' . For the radiation period, 
one has that Xh = — (teq. t) l ' 2 and djj{t) = 2t = H l . We see, for the matter dominated 

era, the causal horizon is twice the Hubble distance, H (sometimes called Hubble 
horizon), and they are equal to each other during the radiation dominated era; therefore, 
one uses them interchangeably. It is clearly seen for both eras that as t — > 0, the Universe 
is causally disconnected, being a(t) > dn(t). But, on the other side, by that time the 
CMBR was already highly isotropic and with a black body spectrum. Then, one has to 
take for granted that the initial conditions for all small horizon volumes were very fine 
tuned to account for the present observed large angle CMBR levels of isotropy, with 
8T jT rs few x 10 5 . This is the horizon problem. 

Event Horizon. The event horizon, d e , determines the region of space which will 
keep in causal contact after some time; that is, it delimits the region from which one can 
ever receive (up to some time t mSiX ) information about events taking place now (at the 
time t): 

ftmux At' 

d e (t)=a(t) — (23) 

Jt ciyt ) 

For a flat model during its matter dominated era (a ~ t 2 ' 3 ), d e — > °° as t max — > °°. 

Luminosity distance. The luminosity distance is the distance measured using the 
energy flux (J^) observed by a light source with absolute luminosity (L s ): 

d 2 L = -^- , (24) 

where & = Lq/S, being Lq the observed luminosity and S = 4n(aofk(x)) 2 is the sphere 
area at z = 0. The luminosity distance becomes 

4=(aofk(x)) 2 r-- (25) 

If we express the energy emitted by a light pulse in a time interval Ati as AE\, the 
absolute luminosity is given by L s = AEi/Ati. Similarly we define Lq = AEo/Ato, where 
A£o is the detected energy in a time AtQ. On the other hand, since the photon energy can 
be expressed in terms of its wavelenght (A), one has that AEi/AEq = Xq/X\ — 1 +z and 
moreover c — 1, being constant implies that X\/At\ = Xo/Ato, from which we finally 
have that 

and the luminosity distance becomes 

d L = a f k (x)(l+z). (27) 



Since fk(x) depends on % and this on the redshift, cf. Eq. (21), thus by measuring 
the luminosity distance, we can determine the expansion rate of the Universe. For the 
ACDM model one finds that [16] 

(i+z) r * (28) 



H « h Jl^h+z')^'- 



The luminosity distance becomes larger when the cosmological constant is present and 
this is what was found to fit better the supernovae la. 

In practice, one uses the relationship of the apparent (m) and absolute (M) magnitude, 
related to the luminosity measured at present and when emitted, respectively to have: 

m-M = 5Lo gK) Q^+25. (29) 

By adjusting best fit curves to their data, two different supernova groups [17] found 
a clear evidence for A in the late 90's. The presence of a cosmological constant makes 
the Universe not only expanding, but accelerating and, in addition, its age is older, not 
conflicting with globular cluster ages. With the course of the years, various supernova 
groups have been getting more confidence that the data is compatible with the presence 
of dark energy, dark matter, and a high value of Hubble parameter. One of the latest data 
release, the Union2 compilation [18], reports that the flat concordance ACDM model 
remains an excellent fit to the data with the best fit constant equation of state parameter 
w = -0.997+0054 for a flat Universe, and w = -1.035+0059 wim curvature. Also, they 
found that Q. m = 0.270 ±0.021 (including baryons and DM) for fixed Q,^ = 0. That is, 
Q A = 0.730 ±0.021. 

Angular diameter distance. The angular diameter distance is given by 

Ax 
d A = Jq, (30) 

where A0 is the angular aperture of an object of size Ax orthogonal to the line of sight 
in the sky. Usually this distance is used in the CMBR anisotropy observations, since the 
source emitting the radiation is on a surface of a sphere of radius % with the observer 
located at the center; it is also used in the determination of the BAO feature, see below. 
Thus, the size Ax at the time t\ in the Friedmann metric, Eq. (19), is given by 

Ax = a(t l )Mx)Ae. (31) 

Thus, the angular diameter distance is 

d A = a{t l )f k { X ) = a i^ ( 32 ) 

1 +z 



and comparing it with Eq. (27) one has 

dh 

which is called duality relationship. Eq. (33) is valid beyond the FRW metric. In fact, it 
is valid for any metric in which the flux is conserved. 



d A = ,-, , x 2 , (33) 



THE PHYSICAL UNIVERSE 

In the following we provide with theoretical tools to understand the physics of the early 
Universe. We treat some micro-physics that rules the interactions of the particles and 
fields. 



Thermodynamics in the early Universe 

In the early Universe one considers a plasma of particles and their antiparticles, as 
originally was done by Gamow [1], who has first considered a physical scenario for 
the hot big bang model for the Universe's beginning. Later on, with the development 
of modern particle physics theories in the 70's it was unavoidable to think about a 
physical scenario which should include the "new" physics for the early Universe. It was 
also realized that the physics described by GR should not be applied beyond Planckian 
initial conditions, because there the quantum corrections to the metric tensor become 
very important, a theory which is still in progress. So the things, one assumes at some 
early time, t^tpi, that the Universe was filled with a plasma of relativistic particles which 
include quarks, leptons, and gauge and Higgs bosons, all in thermal equilibrium at a very 
high temperature, T, with some gauge symmetry dictated by a particle physics theory. 

Theoretically, in order to work in that direction one introduces some thermodynamic 
considerations necessary for the description of the physical content of the Universe, 
which we would like to present here. Assuming an ideal-gas approximation, the number 
density n, of the particles of type /, with a momentum q, is given by a Fermi or Bose 
distribution [19]: 

2dq (34) 



' In 2 J e (E,-H)/T±i ' 

where E t = Jmj + q 2 is the particle energy, jU, is the chemical potential, the sign (+) 

applies for fermions and (— ) for bosons, and g, is the number of spin states. One has 
that gi = 2 for photons, quarks, baryons, electron, muon, tau, and their antiparticles, 
but gj = 1 for neutrinos because they are only left-handed. For the particles existing 
in the early Universe one usually assumes that \X[ = 0: one expects that in any particle 
reaction the /i, are conserved, just as the charge, energy, spin, and lepton and baryon 
number, as well. For photons, which can be created and/or annihilated after some 
particle's collisions, its number density, riy, must not be conserved and its distribution 
with \Xy = 0, E = q = hco, reduces to the Planckian one. For other constituents, in order 
to determine the /i ; , one needs nf, one notes from Eq. (34) that for large /i ; > 0, ni is 
large too. One does not know m, but the WMAP data constrains the baryon density at 
nucleosynthesis such that [20]: 

^ _ ns _ "baryons -"anti-baryons = g 14±0 25 x 1Q -10 (35) 

fly fly 



The smallness of the baryon number density, ng, relative to the photon's, suggests that 
"leptons ma y also be small compared with riy. Therefore, one takes for granted that /i, = 
for all particles. Why the ratio ns/ny is so small, but not zero, is one of the puzzles of 
the standard model of cosmology called baryogenesis, as we will explain below. 

The above approximation allows one to treat the density and pressure of all particles as 
a function of the temperature only. According to the second law of thermodynamics, one 
has [13]: dS(V, T) = j[d(pV) +pdV], where S is the entropy in a volume V ~ a 3 (t) with 
p = p(r), p = p(T) in equilibrium. Furthermore, the following integrability condition 

WW = dvdr * s a ^ so va lid> which turns out to be 

% = *¥■ <36) 

dT T 

On the other hand, the energy conservation law, Eq. (6), leads to 

a3{t) t = Jt [a3{t){p+p)] (37) 

and using Eq. (36), the latter takes the form ^[^^(p + p)] = 0, and using Eq. (36) 
again, the entropy equation can be written as dS(V, T) = jd[{p + p)V] — jj (P + p)dT . 
Last two equations imply that the entropy is a constant of motion: 

a 3 
S=-[p+/>]= const. (38) 

The density and pressure are given by 

p = / Emidq , p= jj^wdq ■ ( 39 ) 

For photons or ultra-relativistic fluids, E = q, these equations become such that p = \ > p, 
and then confirming Eq. (7) with w — 1/3, and after integrating Eq. (36), it comes out 
that 

p = bT\ (40) 

with the constant of integration, b. In a real scenario there are many relativistic particles 
present, each of which contributes like Eq. (40). By including all of them, p = Y*iPi 

and p = Y,iPi over ai l relativistic species, one has that b(T) = §q(Nb + %Nf), which 
depends on the effective relativistic degrees of freedom of bosons (/Vg) and fermions 
(/V>); therefore, this quantity varies with the temperature; different /—species remain 
relativistic until some characteristic temperature T ~ m,-, after that the value Nf j (or 
Nb;) contributes no more to b(T). The factor 7/8 accounts for the different statistics the 
particles have, see Eq. (34). In the standard model of particles physics b ~ 1 for T <^i 1 
MeV and b ~ 35 for T > 300 GeV [19]. Also for relativistic particles, one obtains from 
Eq. (34) that 

n = cT 3 , withc=^h/-(N B +^N F ) . (41) 

% l 4 



where £(3) ~ 1.2 is the Riemann zeta function of 3. Nowadays, riy ps ^f T^S' wnere 
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T = r ° 
^2.75 — 2.15°K- 

From Eq. (38), and using the relativistic equation of state given above (w = 1/3), one 
gets that T ~ l/a(t) and from its solution in Eq. (11) one has, 




i ' 



(42) 



a decreasing temperature behavior as the Universe expands. Then, initially at the big 
bang t = U implies T* = °°, the Universe was not only very dense but also very hot. 

The entropy for an effective relativistic fluid is given by Eq. (38) together with its 
equation of state and Eq. (40), S = i b (a T) 3 = const. Combining this with Eq. (42), 

one can compute the value of Mi to be Mi = (^S) 4 ' 3 /b 1 ' 3 ~ 10 116 , since b ~ 35 and 

the photon entropy So = | b (ao Tq) 3 rs 10 88 for the nowadays evaluated quantities 
oo - > dnito) = 10 28 cm and Ty = 2.7 °K. One defines the entropy per unit volume, 
entropy density, to be s = S/V = i?m(NB + ^Nf)T 3 , then, nowadays s ~ ln y . The 
nucleosynthesis bound on rj, Eq. (35), implies that ng/s ~ 10 -11 . 

Now we consider particles in their non-relativistic limit (m ^> T). From Eq. (34) one 
obtains for both bosons and fermions that 

«^(f) 3/ W. 

The abundance of equilibrium massive particles decreases exponentially once they be- 
come non-relativistic; this situation is referred as in equilibrium annihilation. Their den- 
sity and pressure are given through Eqs. (39) and (43) by p = nm and p = nT <^ p. 
Therefore, the entropy given by Eq. (38) for non-relativistic particles, using last two 
equations, diminishes also exponentially during their in equilibrium annihilation. The 
entropy of these particles is transferred to that of relativistic components by augmenting 
their temperature. Hence, the constant total entropy is essentially the same as the one 
given above, but the /—species contributing to it are just those which are in equilibrium 
and maintain their relativistic behaviour, that is, particles without mass such as photons. 
Having introduced the abundances of the different particle types, we would like to 
comment on the equilibrium conditions for the constituents of the Universe, as it evolves. 
This is especially of importance in order to have an idea whether or not a given /—species 
disappears or decouples from the primordial brew. To see this, let us consider nt when 
the Universe temperature, T, is such that (a) T ^> nij, during the ultra-relativistic stage 
of some particles of type / and (b) T <^ m,-, when the particles / are nonrelativistic, both 
cases first in thermal equilibrium. From Eq. (41) one has for the former case that nt ~ T 3 ; 
the total number of particles, ~ n ; a 3 , remains constant. Whereas for the latter case, from 
Eq. (43), n\ ~ T 3 ' 2 e~ m '' T , i.e., when the Universe temperature goes down below m,, the 
number density of the /—species significantly diminishes; it occurs an in equilibrium 



annihilation. Let us take as example the neutron-proton annihilation, one has 

n mp-m n 1.5xl0 10 °K 

-1 ~ e — r— =e t , (44) 

n p 

which drops with the temperature, from near to 1 at T > 10 12 °K to about 5/6 at 
T « 10 11 °K, and 3/5 at T « 3 x 10 10 °K [21]. If this is forever valid, one ends 
without massive particles, and our Universe should have consisted only of radiative 
components; our own existence prevents that! Therefore, eventually the in equilibrium 
annihilation had to be stopped. The quest is now to freeze out this ratio to be n n /n p ^ 1/6 
(due to neutron decays, until the time when nucleosynthesis begins, n n /n p reduces 
to 1/7) in order to leave the correct number of hadrons for later achieving successful 
nucleosynthesis. The answer comes from comparing the Universe expansion rate, H, 
with particle physics reaction rates, T. Hence, for H <T, the particles interact with 
each other faster than the Universe expansion rate, then equilibrium is established. 
For H > r the particles cease to interact effectively, then thermal equilibrium drops 
out. This is only approximately true; a proper account of that involves a Boltzmann 
equation analysis. For that analysis numerical integration should be carried out in which 
annihilation rates are balanced with inverse processes, see for example [22, 19]. In 
this way, the more interacting the particles are, the longer they remain in equilibrium 
annihilation and, therefore, the lower their number densities are after some time, e.g., 
baryons vanish first, then charged leptons, neutral leptons, etc.; finally, the massless 
photons and neutrinos, whose particle numbers remain constant, as it was mentioned 
above. Note that if interactions of an /—species freeze out when it is still relativistic, 
then its abundance can be significant nowadays. 

It is worth to mention that if the Universe would expand faster, then the temperature 
of decoupling, when H ~ T, would be higher, then the fixed ratio n n /n p must be greater, 
and the 4 He abundance would be higher, thus leading to profound implications in the 
nucleosynthesis of the light elements. Thus, the expansion rate cannot arbitrarily be 
modified during the equilibrium era of some particles. Furthermore, if a particle species 
is still highly relativistic (T 3> m,) or highly non-relativistic (T <C mi), when decoupling 
from primordial plasma occurs, it maintains an equilibrium distribution; the former 
characterized by T r a =const. and the latter by T m a 2 =const., cf. Eq. (47). 

There are also some other examples of decoupling, such as neutrino decoupling: dur- 
ing nucleosynthesis there exist reactions, e.g. vv i — > e + e~, which maintain neutrinos 

efficiently coupled to the original plasma (r > H) until about 1 MeV, since yj ~ (jvj^v) • 
Below 1 MeV reactions are no more efficient and neutrinos decouple and continue evolv- 
ing with a temperature T v ~ 1 /a. Then, at T > m e = 0.5 lMeV the particles in equilibrium 
are photons (with Nb = 2) and electron and positron pairs (with Np — 4) to contribute 

2 

to the entropy with b(T) = |q (11/2). Later, when the temperature drops to T <C m e , the 
reactions are no more efficient (r < H) and after the e ± pair annihilation there are only 

photons in equilibrium with b(T) = §q(2). Since the total entropy, S = ^b(aT) 3 , must 
be conserved, the decrease in b(T) must be balanced with an increase in the radiation 

temperature, then one has that j? = (^) , which should remain so until today, im- 
plying the existence of a cosmic background of neutrinos with a temperature today of 



T Vo = 1.96 °K. This cosmic relic has not been measured yet. 

Another example of that is the gravitation decoupling, which should be also present 
if gravitons were in thermal equilibrium at the Planck time and then decouple. The 

1 1% 

today background of temperature should be characterized at most by r grav . = {jqj) ~ 
0.91 °K. 

For the matter dominated era we have stressed that effectively p = 0; next we will see 
the reason of this. First consider an ideal gas (such as atomic Hydrogen) with mass m, 
then p =nm+ \riT m and p = nT m . From Eq. (37) one obtains, equivalently, that 

^-(pa 3 (t)) = -3pa 2 (t) (45) 

da 

and substituting the above p and p, one has that 

— (nma 3 (t) + -nT m a 3 {t)) = -3nT m a 2 (t) (46) 

where nma 3 (t) is a const. This Eq. yields that 

T m a 2 (t) = const. , (47) 

the matter temperature drops faster than that of radiation as the Universe expands, cf. Eq. 
(42). Now, if one considers both radiation and matter, one has that p = nm + \riT m + bT 4 
and p = nT m + \bT 4 ; the source of Universe's expansion is proportional to p + 3p = 
nm + \riT m + 2bT 4 , the first term dominates the second, precisely because T m decreases 
very rapidly. The third term diminishes as ~ 1 /a 4 , whereas the first as ~ 1 /a 3 , and after 
the time of densities equality, p m = p r , the matter density term is greater than the others, 
that is why one assumes no pressure for that era. 

From now on, when we refer to the temperature, T, it should be related to the ra- 
diation temperature. The detailed description of the Universe thermal evolution for the 
different particle types, depending on their masses, cross-sections, etc., is well described 
in many textbooks, going from the physics known in the early 70's [13] to the late 80's 
[19], and therefore it will not be presented here. However, we notice that as the Uni- 
verse cools down a series of spontaneous symmetry-breaking (SSB) phase transitions 
are expected to occur. The type and/or nature of these transitions depend on the specific 
particle physics theory considered. Among the most popular ones are Grand Unifica- 
tion Theories (GUT's), which bring together all known interactions except for gravity. 
One could also be more modest and just consider the standard model of particle physics 
or some extensions of it. Ultimately, one should settle, in constructing a cosmologi- 
cal theory, up to which energy scale one wants to describe physics. For instance, at a 
temperature between 10 14 GeV to 10 16 GeV the transition of the SI/ (5) GUT should 
took place, if this theory would be valid, in which a Higgs field breaks this symmetry 
to SU(3)c x SU(2)w x U(1)hc> a process through which some bosons acquired their 
masses. Due to the gauge symmetry, there are color (C), weak (W) and hypercharge 
(HC) conservation, as the subindices indicate. Later on, when the Universe evolved to 
around 150 GeV the electroweak phase transition took place in which the standard model 
Higgs field broke the symmetry SU(3) C x SU(2) W x U(l) HC to SU(3) C x U(1)em', 



through this breaking fermions acquired their masses. At this stage, there were only 
color and electromagnetic (EM) charge conservation, due to the gauge symmetry. Af- 
terwards, around a temperature of 175 MeV the Universe should underwent a transition 
associated to the chiral symmetry-breaking and color confinement from which baryons 
and mesons were formed out of quarks. Subsequently, at approximately 10 MeV begun 
the synthesis of light elements (nucleosynthesis), when most of the today observed Hy- 
drogen, Helium, and some other light elements abundances were produced. So far the 
nucleosynthesis represents the earliest scenario tested in the standard model of cosmol- 
ogy. After some thousand years (z ~ 3200 [10]), the Universe is matter dominated, over 
the radiation components. At about 380, 000 years (z ~ 1090 [10]) recombination took 
place, that is, the Hydrogen ions and electrons combined to compose neutral Hydrogen 
atoms, then matter and EM radiation decoupled from each other; at this moment (bary- 
onic) matter structure begun to form. Since that moment the surface of last scattering 
of the CMBR evolved as an imprint of the early Universe. This is the light that Penzias 
and Wilson first measured, and was later measured in more detail by BOOMERANG, 
MAXIMA, COBE, and WMAP, among other probes. PLANCK cosmological data will 
be forthcoming in 2012. 



Inflation: the general idea 

As we mentioned above, the FRW cosmological Eqs. (4)-(6) admit very rapid ex- 
panding solutions for the scale factor. This is achieved when p + 3p, is negative, i.e., 
when the equation of state admits negative pressure such that w < —1/3, to have a > 0. 
For instance, if w = —2/3, one has that a ~ t 1 and p ~ I /a, that is, the source of rapid 
expansion decreases inversely proportional with the expansion. Of special interest is the 
case when w = — 1, p = const., because this guarantees that the expansion rate will not 
diminish. Thus, if p = const, is valid for a period of time, T, the Universe will experience 
an expansion of N = xH foldings, given by a = a*e N , Eq. (12). This is the well known 
de Sitter cosmological solution [23], achieved here only for a T-stage in a FRW model. 

We shall now see how an inflationary stage helps to solve the horizon and flatness 
problems of the old standard cosmology. Firstly consider the particle (causal) horizon, 
given by Eq. (22), during inflation, again with k = 0, one obtains 

d H =H-\e Ht -l) , (48) 

the causal horizon grows exponentially, whereas H remains constant. We compare the 



horizon distance with that of any physical length scale, L{t) = L*— = L*e Ht , to get 

l-e~ m , (49) 



d H _H- l (e Ht -l) >i „ t 



for initial length scales L*<H l . After a few e-fold times the causal horizon is as big 
as any length scale that was initially subhorizon sized. Therefore, if the original patch 
before inflation is causally connected, and presumably in equilibrium, then after inflation 
this region of causality is exponentially bigger than it was, and all the present observed 



(apparent) Universe can stem from it, solving the horizon problem. In fact, if the inflation 
stage is sufficiently large, there can exist nowadays regions which are so distant away 
from each other that they are still not in contact, even though originally they come from 
the same causal patch existing before inflation. These regions will be for a time not in 
contact since we currently are experiencing an accelerated expansion, and then, if this is 
of exponential type, the event horizon is constant and light/information that shall come 
to us will be from only a delimited region H , as we explain below. 

From Eq. (49) one can observe that if the initial physical length scale is greater than 
the Hubble distance, L* > H , then du < L during inflation. Events initially outside 
the Hubble horizon remain acausal. This is better seen by considering the event horizon, 
d e , defined in Eq. (23). This delimits the region of space which will keep in causal 
contact after some time; that is, it delimits the region from which one can ever receive 
(up to some time t mSiX ) information about events taking place now (at the time t). During 
inflation one has that 

rf 8 = r 1 (i-r^ _, ' ff )«/r 1 ) (50) 

which implies that any observer sees only those events that take place within a distance 
< H . In this respect, there is an analogy with black holes, from whose surface 
no information can get away. Here, in an exponential expanding Universe, observers 
encounter themselves in a region which apparently were surrounded by a black hole 
[24, 25], since they receive no information located farther than H . 

The apparent horizon at present stems from a region delimited by the original patch 
d e ~ H , which during inflation remains almost constant and, afterwards, evolves as 
H l ~ t. At the end of inflation a(t) ^> H~^(t)/HZ . Subsequently, the scale factor 

expands only with the power law solution t 1 ' 2 (and later as t 2 ' 3 ), whereas the Hubble 
horizon evolves faster, H ~ t. Then, at some later time the Hubble horizon is as large 
as the scale factor, H l ~ a(t)HQ l . Accordingly, there is a minimal number of efolds 
of inflation, N ~ 60, necessarily to have this equality at present (this number depends 
on the energy scale of inflation, see for instance J. L. Cervantes-Cota in [5]); that is, the 
original patch grown until now is as big as our apparent, Hubble horizon. Hence, some 
time ago, say, at the last scattering surface (photon decoupling) the Universe consisted 
of 10 5 Hubble horizon regions, yet all these regions stem from one original patch of size 
H^ 1 at the start of inflation. 

A typical scale L* < H~ l will increase exponentially its size as L(t) = L* — = L*e N . 
That is, all physical inhomogeneities, anisotropics and/or 'perturbations' of any kind 
(including particles!) will be diluted away from a region d e ~ H~ , and its density 
becomes insignificant, thus solving the monopole (and other relics) problem. 

On the other hand, the flatness problem in the old standard cosmology arises since Q. 
approaches closely to unity as one goes back in time in a way that one has to choose 
very special initial density values, at the Planck time Q. P i — 1 f» ±10~ 59 , for explaining 
our flatness today, i.e., Q.q ps ^(1). Now, imagine the Universe with initial conditions 
such that Q* — 1 ?a k. Now, if the exponential expansion occurs, Q(t = t) evolves to 

AM" 1 = ^ = 34* = *"-"'■ <*) 



If N is sufficiently large, which will be case since typically N > 60, the Universe looks 
after a de Sitter stage like an almost perfect flat model. Therefore, it plays almost no 
role what the initial density was, if the exponential expansion occurs -guaranteed by 
the cosmological no hair- the Universe becomes effectively flat. In this way, instead of 
appealing to very special initial conditions, one starts with an Universe with more normal 
conditions, that is, non fine tuned, which permit the Universe to evolve to an inflationary 
stage, after which it looks like it would had very special conditions, i.e., with Q. as 1 with 
exponential accuracy. 

After inflation the Universe contains a very small particle density and is very cold, 
even as cold as the CMBR is today! The transition to a radiation dominated era with 
sufficient entropy and particle content comes from the 'decaying' or transformation of 
the energy source of inflation, p = V (0), into heat; a process called reheating (RH). 



Reheating and baryogenesis 

At the end of inflation the <j> -field (the inflaton) begins to oscillate around its stable, 
global minimum, say v. Its oscillation frequency is given by the effective mass of the 
Klein Gordon equation in a FRW Universe, + 3H<j) + V'((j)) = 0. After inflation the 
term V'(~ V"ty = M 2 0) is greater than 3H(j), since the Hubble rate evolves hereafter 
always decreasing, H ~ \/t. Thus, the oscillation frequency of the 0-field is simply 
given by the field mass, M. The stored energy of the inflaton field, p# = V + \§ 2 , can 
decay to give rise to quantum particle creation [26]. The reheating models depend on 
the particle physics models but general features of the process have been understood. 
We firstly explain the old scenario called reheating and secondly the preheating that 
seems to be more realistic. In reheating, the state = v is considered as a coherent 
state of scalar particles in rest. Then, this state decays through the ordinary decay of 
the field bosons and the decay rate coincides with the rate of decrease of the energy of 
oscillations. Thus, the decay rate of the boson, T#, introduces a friction term of the type 
Th<J> in the Klein Gordon equation that causes the scalar field to vanish and the reheating 
of the Universe [26]. The transformed energy goes into masses and kinetic energy of 
the new particles produced. Typically, the produced particles (bosons, fermions) have 
smaller masses than the field boson, therefore much of their energy goes into kinetic 
energy, and particles behave as a relativistic fluid. If the decay rate is greater than the 
Hubble rate after inflation, T# > H(t = tf) = Hf, then the reheating process occurs 
within an expansion time, very rapidly. In this case, the reheating temperature is [27] 

Trh « {f) 1 / 4 = (^-j , where b = ^(N B + 7 E N F ), N B (N F ) stands for boson 

(fermion) degrees of freedom. For high temperatures, T > 300 GeV, one has that b ~ 35, 
see discussion after Eq. (40). The subindex / means to be evaluated at the end of 
inflation. There, the kinetic and potential energies are of the same order of magnitude. 
Then, for a self-interaction potential with V « V(0) = ^Av 4 , one has that 



iRH 



Av^ 1/4 
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(52) 



One the other hand, if the decay process is rather slow, r# < Hf, then the field 
continues to oscillate coherently until t = T^ 1 . During this time the solution of the Klein 
Gordon equation is <j) ~ j -cosMt, H = ^. The coherent oscillations behave as non- 

relativistic matter fluid (w = 0), i.e., a ~ ? 2 ' 3 [28]. If field bosons do not completely 
decay, the oscillations represent a "sea" of cold bosons with M ^> T . They can account 
for the cold dark matter, but some degree of fine tuning is necessary [29]. 

Without partial or total decaying of field oscillations the Universe remains cold 
and devoid of fermions and (other) bosons. Therefore, let us suppose that indeed re- 
heating took place, but now with T^^Hf, then the reheated temperature is Trh rs 
(\/Mv/(\^Sb) ^Tn/Hf, a factor ^Tn/Hf smaller than the efficient reheating case, 
Eq. (52). The reheating process occurs normally within one or few Hubble times. Then, 
the scale factor does not increase significantly during it. 

The reheating scenario presented above is based on the original theory developed 
in the context of the new inflationary scenario, however, it is also applicable to other 
models. In the course of the time important steps to consolidate the theory were made, 
see for example Ref. [30]. But qualitative new ideas were introduced in Refs. [31]. 
Accordingly, the process of reheating should consist of three different stages. At the first 
phase, the (j) -field decays into massive bosons (fermions) due to a parametric resonance 
given through a Mathieu equation that determines the regions of stability and instability 
(particle production) of the quantum fluctuations of the created particles. These can be 0- 
particles or other bosons (fermions) coupled to the 0-field. This process is very efficient, 
even explosive, and much bosons can be created in this stage. Note that the original 
theory is based upon the decay of the -particles, whereas in the present theory the <j>- 
field decays into ^-particles, and perhaps others, and only after this process the decay 
of these particles proceeds. Then, to distinguish this explosive process from the normal 
stage of particle decay, the authors of Ref . [31] called it preheating. Bosons produced 
at this stage are far away from thermal equilibrium and have very big occupational 
numbers. The second stage of this scenario describes the decay of the already produced 
particles. This phase is described as in the original reheating theory. Then, the methods 
developed for the original theory are now applied to the product particles, but not itself 
to the decay of the <j> -field. The third stage is the thermalization by which the system 
reaches equilibrium [32]. 

The process of reheating is very complex and depends fine on the particle physics 
theory one has in turn. As a matter of fact, one expects a reheat temperature Trh > 
few MeV to be able to attain nucleosynthesis. A second, and more restrictive constraint 
comes from baryogenesis. One can see this by noting that the number density of any 
conserved quantity before reheating divided by the entropy density - becomes after 
reheating insignificant because of the huge entropy produced. One gets "| , = e~ 3N ^,\ ri , 
where ri and rf denote the initial and final state of reheating, respectively. In this 
way, any baryon asymmetry initially present will be brought to unmeasurable values. 
Therefore, after reheating the baryon asymmetry must be created. Note also that any 
unwanted relic (x), accounted through — , will essentially disappear after reheating. The 
correct baryon-antibaryon balance at nucleosynthesis is ng/s ~ 10~ u , in consistency 
with Eq. (35). There are some attempts to achieve baryogenesis at low energy scales, as 
low as few GeV or TeV [33]. Recent attempts to solve this problem seek to yield a prior 



a lepton asymmetry, leptogenesis, generated in the decays of a heavy sterile neutrino 
[34], to later end with baryogenesis. 

THE PERTURBED UNIVERSE 

In the previous sections we have outlined how the evolution of a homogeneous Universe 
can be described by means of few equations and simple concepts such as the ideal perfect 
fluids. The next step is introducing in this scenario small inhomogeneities that can be 
treated as first order perturbations to those equations, the goal being the description of 
the structures we see today in the Universe. This perturbative approach is sufficient to 
accurately describe the small temperature anisotropics (AT /T ~ 10~ 5 ) observed in the 
CMBR today, but can describe the distribution of matter today only at those scales that 
are still in the linear regime. At the present epoch, scales smaller than ~ 30 Mpc hr x 
[35] have already entered the non linear-regime (Ap/p >> 1) due to the fact that matter 
tends to cluster under the effect of gravity. These scales can therefore only be described 
by means of numerical or semi-numerical approaches [36]. 

The approach is quite straightforward but involves a differential equation for the 
density perturbation of each individual constituent: scalar fields in inflation, or baryons, 
radiation, neutrinos, DM, and DE (usually treated as cosmological constant) in later 
times, and in general needs to be solved numerically. In the context of geometry metric 

and/or GR the metric is treated as the general expansion term g^{, plus perturbation h^ v : 

g{iv = gfl + h iiV , (53) 

with h^v << g^y where (°) indicates the unperturbed homogeneous quantities. 

Inhomogeneities in the distribution of the components of the Universe are a source of 
scalar perturbations of the metric. Nevertheless vector or tensor perturbations can modify 
the metric as well. The standard cosmological model do not predict vector perturbations, 
that would introduce off-diagonal terms in the metric tensor. These perturbations would 
produce vortex motions in the primordial plasma that are expected to rapidly decay. 
Models with topological defects or inhomogeneous primordial magnetic fields instead 
predict a consistent fraction of vector perturbations [37, 38, 39]. 

On the other hand, the standard cosmological model predicts the production of gravi- 
tational waves during the epoch of inflation, when the Universe expanded exponentially, 
as we will see below. Gravitational waves induce tensor perturbations hL v on the metric 
of the type: 

/ \ 

h+ h x 

h x -h+ 
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h T =a 2 



where h + and h x are the polarization directions of the gravitational wave. This tensor 
is traceless, symmetric and divergentless, i.e. it perturbs the time space orthogonally 
to the propagation direction of the wave. The amplitude of these tensor perturbations 
are expected to be small compared to scalar ones, and therefore negligible in first 



approximation as far as we are interested in studying the perturbations of the metric 
tensor. Nevertheless these waves are expected to leave an imprint in the polarization 
of the CMBR, and their eventual detection would unveil an extremely rich source of 
information about an epoch of the Universe that is very hardly observable otherwise. 

It is important to underline that choosing to model the metric perturbations corre- 
sponds to choosing a gauge, i.e. a specific coordinate system in which the metric tensor 
is represented. Changing the coordinate system of course do not change the physics, but 
can remarkably vary the difficulty of the calculations and the understanding of the phys- 
ical meaning of the different quantities. To solve the perturbed equations one chooses 
convenient gauges for the different expansion epochs and depending on whether the 
formalism is theoretical or numerical, as we will see below. 

The presence of weak inhomogeneous gravitational fields introduce small pertur- 
bations in the metric tensor that can be modeled by introducing two scalar functions 
<J>(x, T])) and ^(x , rj) [40] in the Robertson- Walker metric as: 

ds 2 = a 2 (ri) [-[1 +2<S>(x,ri)}dri 2 + [1 +2 x ¥(x,ri)}dx i dx i ] , (54) 

where therefore the perturbed part of the metric tensor is: 

h o(x,ri) = -20>(x,T]), h 0i (x,ri) = 0, hij(x,r\) = a 2 8 ij (2 x ¥(x,ri)). (55) 

This metric is just a generalization of the well known metric for a weak gravitational 
field usually presented in text books (e.g. Chapt. 18 in Misner [7]) for the case of 
a static Universe (a(rj) = 1). The function 4> describes Newton's gravitational field, 
while *F is the perturbation of the space curvature. The above gauge is the Newtonian 
conformal gauge, which has the advantage of having a diagonal metric tensor g^ v in 
which the coordinates are totally fixed with no residual gauge modes and therefore with 
a straightforward interpretation of the functions introduced [40, 41, 42]. An example of 
an alternative gauge particularly popular in literature is the synchronous gauge, generally 
defined as: 

ds 2 = a 2 (r])[-dri 2 + (8ij + hij)dx i dx j ], (56) 

which is especially used in codes computing the anisotropics and inhomogeneities in the 
Universe, as better behaved numerically by choosing that observers fall freely without 
changing their spatial coordinates. A further analysis is found in e.g. [43]. 



Perturbations during inflation 

The primeval fluctuations are thought to be present at the very beginning of time, at 
the inflationary epoch. The generation of perturbations are produced by quantum fluc- 
tuations of the -field during the accelerated stage, for a review see [40, 5, 43]. These 
fluctuations are usually studied in the comoving gauge in which the scalar field is equal 
to its perturbed value at any given time during inflation and, therefore, the perturbation 
information resides in the metric components. The perturbations cross outside the event 
horizon during inflation and re-enter into the horizon much later, at the radiation and 



matter dominated epochs, to yield an almost scale invariant density perturbation spec- 
trum (Harrison-Zel'dovich, n$ = 1), as the required for structure formation. 

We introduce this topic by noting that the event horizon during a de Sitter stage is 
d e ~ H , cf. Eq. (50). This means that microphysics can only operate coherently within 
distances at most as big as the Hubble horizon, H . Recall that the causal horizon, dn, 
expands exponentially and it is very large compared to the almost constant H during 
inflation, see Eq. (48). Hence, during the de Sitter stage the generation of perturbations, 
which is a causal microphysical process, is localized in regions of the order of H . 

It was shown that the amplitude of inhomogeneities produced corresponds to the 
Hawking temperature in the de Sitter space, 7# = H /{2k). In turn, this means that 
perturbations with a fixed physical wavelength of size H are produced throughout 
the inflationary era. Accordingly, a physical scale associated to a quantum fluctuation, 
■^phys = ka{t), expands exponentially and once it leaves the event horizon, it behaves as 
a metric perturbation; its description is then classical, general relativistic. If inflation 
lasts for enough time, the physical scale can grow as much as a galaxy or horizon 
sized perturbation. The field fluctuation expands always with the scale factor and after 
inflation, it evolves according to t n in = 1/2 radiation or n = 2/3 matter). On the other 
hand, the Hubble horizon evolves after inflation as H ~ t. This means, it will come a 
time at which field fluctuations cross inside the Hubble horizon and re-enters as density 
fluctuations. Thus, inflation produces a gross spectrum of perturbations, the largest scale 
originated at the start of inflation with a size Hf , and the smallest with HJ X at the end 
of inflation. The power spectra for scalar (5) and tensor (T) perturbations are given by: 



ft(i)w (A)L- witzk 



(57) 



k—aH 



where <j) c is the classical field velocity. The equations are evaluated at the horizon 
crossing (k = aH) during inflation. Each of the k— modes generate also an anisotropy 
pattern in the CMBR that was measured for scalar perturbations by the COBE [4] 
and later probes. The PLANCK satellite may have the chance to detect the ratio of 
tensor to scalar amplitudes r = Cf /Cf < 0.36 [10], since the tensor modes modulate 
CMBR photons coming from last scattering. Associated to these perturbations one has 
the spectral indices, n$ ~ 0.96 [10] and nj, and their runnings, dng/dlnk ~ —0.034 [10]. 
These density and other metric perturbations are small, but we discuss in the next 
section how to include them so that the information contained can be recognized and 
exploited. 



Perturbations inside the horizon 

In the early Universe, baryons were tightly coupled to photons in an expanding back- 
ground. Baryonic matter and dark matter potential wells provoked the local collapse of 
density fluctuations up certain point, at which the radiation pressure was big enough to 
pull out the matter apart, and smoothing the potential wells. These oscillations of the 
plasma can be thought of as acoustic waves. As we know any wave can be decomposed 
into a sum of modes with different wave numbers, k = 2it/X. Since these modes are 



in the sky, their wavelengths are measured as angles rather than as distances. Accord- 
ingly, instead of decomposing the wave in a Fourier series, what is normally done is to 
decompose the wave in terms of spherical harmonics, Y[ m (p). The angular power spec- 
trum can be expanded in Legendre polynomials, since there is no preferred direction in 
the Universe and that only angular separation 6 is relevant. A mode / plays the same 
role of the wavenumber k, thus / ~ 1/0. Ultimately, we are interested in the tempera- 
ture fluctuations that are analyzed experimentally in pairs of directions h and h' where 
cos (6) = h-h' . We then average these fluctuations obtaining the multipole expansions: 

AT °° l (21 + 1) 

y = n a lm (x,r,)Y lm (p), Ps (e) = J £ K -^— L C,P l (cose), (58) 

l=lm=—l 

where Ps(6) is the angular power spectrum, Pi are the Legendre polynomials and Q is 
estimated as the average over m of a/ m . All this information can be used to determine 
the cosmological parameters Q ; . We will not discuss detailed calculations nor the curve 
that must be adjusted to obtain the best fit values for such parameters. The peak of the 
fundamental mode appears at approximately 

200 

'-7E- (59) 



BOOMERANG [44] and MAXIMA [45] were two balloon-borne experiments de- 
signed to measure the anisotropics at smaller scales than the horizon at decoupling 
($hor-dec ~ 1°), hence measuring the acoustic features of the CMBR. The sensitivity 
of the instruments allowed a measurement of the temperature fluctuations of the CMBR 
over a broad range of angular scales. BOOMERANG found a value of / = 197 ± 6 and 
MAXIMA- 1 found a value of / ~ 220. This implies that the cosmological density pa- 
rameter Q. rs 1, see Eq. (9), implying that the Universe is practically flat, Q.^ ps 0. These 
two experiments provided with the first strong evidence for a flat Universe from observa- 
tions. Happily, this result was expected from Inflation. These results were confirmed by 
the Wilkinson Microwave Anisotropy Probe (WMAP) in a series of data releases in the 
last decade, as well as by other cosmological probes: the Universe is flat or pretty close 
to be flat. The problem in the exact determination of the curvature is because the CMBR 
anisotropics show strong degeneracies among the cosmological parameters [46]. How- 
ever, the satellite PLANCK will offer results of the density parameters with uncertainties 
less than a percentage level. 

Since baryons and photons were in thermal equilibrium until recombination (also 
called last scattering), the acoustic oscillations were imprinted too in the matter per- 
turbations, as they were in the CMBR anisotropics. These are known as baryon acoustic 
oscillations (BAO). The sound horizon at the moment when the baryons decoupled from 
photons plays a crucial role in determination of the position of the baryon acoustic peaks. 
This time is known as drag epoch which happens at Zd = <zo/°d ~~ 1- The sound horizon 
at that time is defined in terms of the effective speed of sound of the baryon-photon 

V 



plasma, c 2 s = 8py/(8py + 8pb) 



fVd 1 r a d da 

r s (z d )= dric s {ri) = - /rrr „ //10 ? • (60) 

Jo 3 Jo a z H(a)y/l + (3Q. h /4Q.y)a 



Note that drag epoch does not coincide with last scattering. In most scenarios z<\ < zi s 
[47]. The redshift at the drag epoch can be computed with a fitting formula that is a 

function of Q. m h 2 and Q. h 'h 2 [48]. The WMAP -5 year team computed these quantities 
for the ACDM model obtaining z d = 1020.5 ± 1 .6 and r s (z d ) = 153, 3 ± 2.0 Mpc [49]. 
What one measures is the angular position and the redshift [50, 6]: 

S (Z) = J^Mj-y (61) 

{l+z)d A {z) 

8z s (z) = r s {z d )H(z), (62) 

where d A (z) is the proper (not comoving) angular diameter distance, Eq. (32), and H(z) 
by Eq. (16). The angle 9 s (z) corresponds to the direction orthogonal to the line-of- sight, 
whereas 8z s (z) measures the fluctuations along the line-of-sight. Observations of these 
quantities are encouraging to determine both d A (z) an d H(z). However, from the current 
BAO data is not simple to independently measure these quantities. This will certainly 
happen in forthcoming surveys [51]. Therefore, it is convenient to combine the two 
orthogonal dimensions to the line-of-sight with the dimension along the line-of-sight to 
define [52]: 

/ 7 \ 1/3 

D v (z)^^l+z) 2 d A (z) 2 1 ^ ) j , (63) 

where the quantity Dm = d A /a = (1 +z)d A (z) is the comoving angular diameter distance. 
One also defines the BAO distance 

rBAo(z) = r s (z d )/D v (z). (64) 

The BAO signal has been measured in large samples of luminous red galaxies from the 
SDSS [52]. There is a clear evidence (3.4a) for the acoustic peak at lOO/z -1 Mpc scale. 
Moreover, the scale and amplitude of this peak are in good agreement with the prediction 
of the ACDM given the WMAP data. One finds that D v (z = 0.35) = 1370 ± 64 Mpc, 
and more recently new determinations of the BAO signal has been published [53] in 
which s (z = 0.55) = 3.90° ±0.38°, and w = -1.03 ±0.16 for the equation of state 
parameter of the dark energy, or Q.m = 0.26 ± 0.04 for the matter density, when the 
other parameters are fixed. 

Measuring the BAO feature in the matter distribution at different redshifts will help 
to break the degeneracy that exists in the determination of the cosmological parameters. 
And by combining line-of-sight with angular determinations of the BAO feature one will 
constrain even more the parameter space. Further, a complete combination of BAO, the 
full matter power spectrum, Supernovae la, and CMBR data shall certainly envisage the 
true nature of the mysterious, dark Universe. 

Finally, we show recent plots of the CMBR power spectrum from the WMAP 7- 
year data and from the South Pole Telescope, from Ref. [54], and the mass variance 
(Am/M = ^P(k)k 3 /(2K 2 ) ) of the reconstructed matter power spectrum from the Ata- 
cama Cosmology Telescope and other observations, from Ref. [55]. The left panel shows 
the extraordinary fit of the ACDM model and the importance of foregrounds for large 
/-modes. The right panel shows the variance decrease as the mass increases, covering 
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FIGURE 1. Left Panel: Recent CMBR angular power spectrum from the WMAP 7-year data and from 
the South Pole Telescope observations show as band-averaged powers along with the best fit ACDM model 
(CMB - dashed line) and (CMB + foreground solid line), taken from Ref. [54]. The detailed location, 
amplitude, and shape of the peaks (bumps) provide information on the contents of the Universe and the 
conditions at that early epoch plus secondary effects. Right panel: The mass variance of the reconstructed 
matter power spectrum from the Atacama Cosmology Telescope; large masses correspond to large scales 
and hence small values of k, taken from Ref. [55]. The BAO are barely visible in the detail of this arching 
spectrum driven by the damping, mostly at low mass, of the original perturbations during the oscillations. 



ten orders of magnitude in the range of masses. We also notice the effect of BAO at 
intermediate scales and damping on the essentially scale invariant perturbations that one 
anticipates from inflation. These observations fit remarkably well to the ACDM model. 
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